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By Melvin F. Mar key 


SUMMARY 


A theory is derived for determining the loads and motions of a 
deeply immersed prismatic body. The method makes use of a two-dimensional 
water -mass variation and an aspect-ratio correction for three-dimensional 
flow. The equations of motion are generalized by using a mean value of 
the aspect-ratio correction and by assuming a variation of the two- 
dimensional water mass for the deeply immersed body. These equations 
lead to impact coefficients that depend on an approach parameter which, 
in turn, depends upon the initial trim and flight-path angles. 

Comparison of experiment with theory is shown at maximum load and 
maximum penetration for the flat-bottom (0° dead-rise angle) model with 
beam- loading coefficients from 36-5 to 133*7 over a wide range of initial 
conditions. A dead-rise angle correction is applied and maximum- load 
data are compared with theory for the case of a model with 30 ° dead-rise 
angle and beam-loading coefficients from 208 to 530 . 


INTRODUCTION 


Several theories exist for the prediction of hydrodynamic- impact 
loads. A survey of the earlier theories which deal primarily with verti- 
cal impacts is presented in reference 1, together with an extension to 
include the effect of velocity parallel to the keel. In reference 2 it 
is shown that this modified theory can be expressed in generalized terms. 
References 1 and 2 deal mainly with the impacts of lightly loaded bodies 
having little or no chine immersion. In reference 3 the theory is 
extended to include the prismatic body with immersed chines, a step-by- 
step procedure being utilized for the determination of the important 
impact parameters. It is the purpose of the present investigation to 
provide a generalized impact theory for the case of deeply immersed 
chines. The formulation of the theory is based on the assumption that 
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the water-mass variation throughout the impact is the same as that for the 
case of deeply immersed chines. Of course, this assumption is not valid 
when the nonimmersed portion of the impact is an important part of the 
total impact, and hence a deeply immersed impact, as far as this report 
is concerned, is one in which the nonimmersed portion plays a relatively 
minor role. In actuality, this condition is approached with the impact 
of a heavily loaded configuration, an important case since the present 
trend is toward such design. 

In the present analysis the equations of motion are written in 
general terms leading to expressions that are applicable to either 
immersed or nonimmersed cases. The expressions are applied to the deeply 
immersed case by the substitution of an assumed two-dimensional water- 
mass variation, resulting in relations that allow such parameters as 
acceleration, penetration, and time to be written in a coefficient form 
that is dependent only upon the velocity ratio and an approach parameter. 
This generalization reduces considerably the number of solutions necessary 
to cover the entire range of impact conditions and also simplifies the 
presentation of experimental results. 

Theoretical and experimental curves of the acceleration and penetra- 
tion coefficients are presented at maximum load and maximum penetration 
for the case of 0° dead-rise angle over a wide range of initial condi- 
tions for bodies having beam- loading coefficients ranging from 36*5 
to 133.7. For bodies having a 30° dead-rise angle and beam-loading coef- 
ficients from 208 to 530, a correcting factor is applied and experimental 
maximum- load data are compared with the theory. Also for the case of 
0° dead-rise angle with beam- loading coefficients of 18.8 and 4.4, cor- 
rection factors are applied to experimental results and these results 
are compared with the present theory for the maximum acceleration 
coefficient. 


SYMBOIS 


b beam 

F hydrodynamic force 

g acceleration due to gravity 

G constant (eq. (l8)) 

l wetted length along keel 

m, two-dimensional water mass in transverse plane 
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s distance along keel from foremost immersed station to flow 

plane 

s velocity along keel 

t time after water contact 

V resultant velocity 

W weight 

x distance parallel to undisturbed water surface, positive in 

direction of body motion 

x horizontal velocity 

x horizontal acceleration 

z immersion of keel at step, measured normal to undisturbed 

water surface, positive downward 

z vertical velocity 

z vertical acceleration 

P angle of dead rise 

7 flight -path angle 

£ distance from undisturbed water surface to keel in any given 

flow plane, measured normal to keel, positive downward 

£ velocity normal to keel 

£ acceleration normal to keel 

A ratio of length of keel below undisturbed water surface to 

the beam 

p mass density 

t trim angle 

cp(A) aspect-ratio correction 

cp^Ai^ mean aspect-ratio correction 
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Subscripts: 

o at water contact 

s at step 

x horizontal direction 

z vertical direction 

£ normal to keel 

max maximum 

Dimensionless parameters: 

beam- loading coefficient, 

sin t 


W 


Pgb' 


approach parameter, ~~ y ■ cos^t + 7 ^ 


acceleration coefficient. 


zb ,£W sin t cos t 


. 1/2 


Ct 


z 0 2 ( G 9(^p) 

l/2 

z ( G ©(An ) pgb^ \ 

penetration coefficient, —I — j 

b \£W sin t cos t / 

ti Q f G 9(^3.) Pgb' 


k3 W 2 


time coefficient, 




b y^W sin t cos t ) 
vertical- velocity ratio 


THEORY 


General 

The basic theory follows closely that of references 1 to 3- It is 
based on the idea that the flow about a slender immersing body occurs In 
transverse planes oriented perpendicular tc the keel. The motion of the 
fluid in each plane is treated as a two-dimensional phenomenon and an 
aspect-ratio correction is applied for the three-dimensional case. The 
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effects of buoyancy, viscosity, and trim changes are neglected and the 
reaction in a given flow plane ds (fig. l) is written as 


dF = 


cLF = 


ikM * 3 


i ^ + ^ ds 


(1) 

( 2 ) 


By integrating equation (2) over the wetted length l the expression 
for the hydrodynamic force is given as 


F = 




(3) 


By using Newton's second law and an aspect-ratio correction qp(A) for 
three-dimensional flow, equation ( 3 ) may be written 



Since 

H Sg _ dry « 
St S£ St d£ 


( 4 ) 


( 5 ) 


and 


d£ = tan t ds 


( 6 ) 


the first term on the right-hand side of equation (4) may be integrated 
and the results rearranged to give 


1 + 


g cp(A) 


W 


r 

d 0 


1 % ds 


g qp(A) 

W tan t 



(7) 


where g is the two-dimensional water mass evaluated at the step. 
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To obtain the equation of motion ( 7 ) in terms of the vertical com- 
ponents, the following substitutions are made: 


i = 


COS T 


and 


i - 


s tan t 


COS T 

Equation ( 7 ) may then be written 


^,5 

sin r 


(z - s sin t)' 


_w__ + r 

<p(?0 J 0 


"V 


ds 


(8) 

( 9 ) 


(10) 


Multiplying equation (l0) by dz, substituting z dz - z dz, and inte- 
grating gives 


L 


z dz 


z D (z - s sin t)‘ 


f 

Jo 


m v, s 

z slrTT dz 


W 


g <p(A) 


+ J "V ds 


( 11 ) 


or 


+ K 


log 


1 + K 


+ K 


1 + K 


r 


^,3 

sin t 


dz 


W 


g 


^*T + 


J 

Jo 


"V ds 


( 12 ) 


where k = - 


s sin t sin t 


z sin 7 n 

and flight-path angles is shown on figure 
in that it may be applied to impacts with 
immersed upon proper substitution of the 
variation and the aspect-ratio correction 
mass variation is considered proportional 
and the aspect-ratio correction is taken 


cos (t + 7 ) and its variation with trim 


2. Equation (12) is general 
chines either immersed or non- 
swo— dimensional water -mass 
. If, for instance, the water - 
to the penetration squared 
as a constant, the equation 
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can be integrated directly, leading to the equations of reference 2. 

In reference 3 an equation similar to equation (12) is derived and applied 

to the immersed- chine case. In the application the term / nu ds is 

J o 

eventually neglected and the equations of motion are solved by a step- 
by-step procedure involving the Pabst aspect-ratio correction and various 
water-mass variations throughout the impact. 

In equation (l2), if it is assumed that cp(A) can be approximated 
by a mean value and taken as a constant with respect to the integration, 
the expression may be integrated directly. The right-hand side of 

equation (12) is of the — form since 

u 


d / W 

dz^g cp(A x ) 



d l 

dz Vs 


m w, s 
sin t 


where cp(Aq) is a representative mean value for the aspect-ratio 
correction. Carrying out the integration of equation (12) leads to 



Evaluation of the right-hand side of equation (13) gives 


+ K 


log 


1 + K 1 + K 


= - log 


\ + g 'P(h) 

r dz] 

W 

Jo sin T 


(14) 


which can be written 


K 

1+k 


g <p(M 

W 


r z V£. dz = l±jn 

J 0 sin T z_ + K 


( 15 ) 


Equation ( 15 ) relates the penetration to the velocity ratio. 
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An expression relating the acceleration, velocity, and penetration 

p 2 

can be obtained by dividing equation (10) by z Q and solving for 

Zq 


The result is 


+ K 


2 g < p(? i) m w J 


W sin t 


z 

. 2 
z o 


g ^p( A i) r z dz 

W J n sin t 


1 + 


(16) 


If the mean value of the aspect— ratio correction can be successfully 
approximated and expressions for the two-dsmensional water mass are avail- 
able, equations (15) and (l6) can be applied in general to get the loads 
and motions of the hydrodynamic impact. It is interesting to note that 
cp(A) is taken as a constant for the nonimmersed-chine case of reference 2 
and it also approaches the constant 1 for the deeply- immersed- chine case 
to be treated. 


Deeply Immersed Chines 


In the application of equations (15) find (l 6) to deeply immersed 
chines, it is necessary to select a variati.on of the two-dimensional 
water mass m w s * In reference 3 ^ha variation suggested is 


mw, s pb 


£(f(0) tan p) : 
o 



tan p 


0 > 



(17) 


where t' represents the normal penetration at the step, including the 
s 

effect of a water rise, and B is a theoretical constant which varies 
with dead-rise angle sometimes called the Bobyleff or Kirchhoff 
coefficient . 

In the present investigation, as an approximation the nonimmersed- 
chine portion of the impact is disregarded and the expression for the 
two-dimensional water mass is written simp.Ly 


,2 5 S h 2 

“w, s = P b G = P b G 


b COS T 


(18) 


It may be noted that this expression is based upon the penetration 
measured from the undisturbed water surfac b and thus neglects the effect 
of water rise. 
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Substituting equation (l8) into equation (15) and integrating gives 

k _ k 

1+K 

Gcp^pgbz 2 1+k 

= e 

2W sin t cos t 7 

— — + K 
Z 0 

Similarly, equation (l6) becomes 



_ 2 /f_ + K \ 2 G 9 (A 1 )pgbz 

z _ \ z o / 2W sin t cos t 

Z Q 2 G cp(A,)pgbz 2 

1 + 

2W sin t cos t 


(20) 


Generalized parameters .- Equations (19) and (20) are similar in form 
to those derived in reference 2 for the nonimmersed-chine case. As in 
reference 2, the equations are written in terms of generalized parameters. 
If, for instance, 


G 9 (A-jJ pgb 


3 v ^ 


c , « - 

Q b\2W sin r cos t. 


( 21 ) 


then from equation (19) 


1/2 


1+K 


- + K 


C d “ 


1+K 


+ K 


and if an acceleration coefficient is defined as 




zb /2W sin t cos t 


,!/2 


G cp(A 1 )pgb' 


then equation (20) may be written 


( 22 ) 


(23) 


C a = 2 f- + k 


1 + C. 


(24) 
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or, in terms of velocity ratio and k , 


C a = 


fz_ 

Uo 


— + K -T-+K 


1+K I 


+ K Z 


1+K 


L + K 



1/2 


( 25 ) 


A time coefficient may also be defined as 

Ct = 


tS o| 

( G cp, 

[A x ) Pgb 9 1 

1 “J 

r 4 dC d 

b ‘ 

l2W sin t cos t J 

0 Z 


( 26 ) 


or 



( 27 ) 


Equations (21) to (27) are the generalized equations of motion for the 
deeply immersed impact. Writing the generalized coefficients in terms 

of the approach parameter k and the vertical-velocity ratio , as 

*o 

indicated, reduces the number of parameters necessary to describe an 
impact . 

Maximum acceleration .- To obtain the ccnditions at maximum accelera- 
tion,^ ” equatTc^ written as 

(l + Kz^)z + 2Kz^z + *z Q )^ = 0 (28) 


where 

G C P(A 1 ) Pgb 

K = ■ — * ' ■ — 

2W sin t cos t 

Equation ( 28 ) can then be differentiated to give 

•zG. + Kz 2 ) + zz( 2 Kz) + 4 -Kzz (z + kz q ) + 2 Kz(z + kz q ) 2 = 0 


( 29 ) 
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At maximum acceleration the differential of the acceleration with respect 
to time, z , is equal to zero; therefore, 


or 


2Kzzz + 4 Kzz (z + kz 0 ) + 2Kz(z + kz q ) - 0 


z 

. 2 



(30) 


(31) 


Equation (31) gives a relationship for the velocity, penetration, and 
acceleration at maximum acceleration. If equation (31) is written in 
coefficient form and combined with equation (24) to eliminate the accelera 
tion coefficient, the following expression is obtained for the penetra- 
tion coefficient at maximum acceleration: 


C 


d 



(32) 


In a similar fashion equation (32) can be used with equation (31) to give 
for the acceleration coefficient: 



(33) 


Equations (32) and (33) can be combined to give the 
the penetration and acceleration coefficients as 


C 


a, max 



relationship between 


( 3*0 



12 


The relation between k and — at maximum acceleration is obtainable 

*o 

from equations (32) and (22). Combining these two equations to eliminate 
leads to 


log 




(35) 


In order to obtain the various generalized coefficients as a function 
of k at maximum acceleration, it is necessary to make a trial-and-error 
solution of equation (35) and then eliminate the velocity ratio from the 
expressions for the generalized coefficients. However, this need be done 
only once . 


Maximum penetration . - The expressions that apply at maximum penetra- 
tion are easily obtained since at this time the vertical velocity ratio 
is equal to zero. Substituting this value Into equation (22) gives for 
the penetration coefficient: 


. 1/2 


'd,max 


1+K 1+K -1 

e - 1, 


acceleration 

coefficient: 



“ / 



i / 

i \ 

C a = 2K 2 

* . e 1+K ll 

- K e 1+K 

1+K \ 

1+K / 


1/2 


(36) 


(37) 


Limiting conditions .- Since the approach parameter k may range 
between 0 and <*, it is of interest to investigate the values of the 
coefficients at these two conditions. In tie former case the flight path 
at initial contact is normal to the keel ani the momentum of the body is 
absorbed entirely by the flow planes normal to it. For this condition 
equation (22) for the penetration coefficient reduces to 



( 38 ) 
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From equation (25) the expression for the acceleration is found to be 


c a . - ^' 1/2 


, \ 


:) 


(39) 


Combining equations ( 38 ) and (39) gives the relationship between the 
acceleration and penetration coefficients: 




2C. 


m 


1 + c. 


For the condition of k = 0 a time coefficient can also be found directly 

by solving equation ( 38 ) for — 9 substituting into equation ( 26 ), and 

*o 

integrating. This gives 


c* 1 + 4 c. 


(41) 


or, in terms of the vertical velocity, 

1/2 


c t = 


a - —\ 


z__ 

7> r 


-A 


1 + 


3 — 


°J 


(* 2 ) 


The equations for the various coefficients at max imum acceleration 
can also be determined for the case of k = 0. From equation (32) it is 
found that 


= 0.2 


This value represents the ratio of the virtual mass to the total mass of 
the body. For the nonimmersed-chine case of reference 2, a value of 2/7 
or 0.286 was obtained for this ratio. Using the value of = o.2 in 

conjunction with equation ( 38 ) gives the value of -5— at maximum acceler- 


L. = 1 = 0.834 

Zo 6 


ation: 
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Reference 2 gives a value of 7/9 = 0-778 far the nonimmersed-chine case. 
The other coefficients may be evaluated as follows: 

C d = \| 0.2 = 0.448 

C a = 2(0.448) _ 0 . 519 

(l + 0.2) 5 

C t = (0.2) 1/2 ^1 + Qy'j = 0.478 

From equation (36) it can be seen that the maximum penetration 
corresponding to the case of k = 0 is ne^er reached as a consequence of 
the neglect of the buoyant forces in the tneory. This may be noted also 
from equation ( 38 ); since letting the vertical-velocity ratio approach 
zero gives a value of the penetration coefficient approaching 00 . 

As k approaches infinity the other end condition is approached, 
that of pure planing. In this instance, the coefficients approach the 
following values: 


lim Ca = 0 

K — > 00 


lim i_ = 1 

* ->°° ; 
z o 


lim C a = 00 

K — > 00 


lim C t = 0 

K — > co 


These results are for the case in which the wing lift is equal to the 
weight . 

For the planing condition with partial wing lift the relationship 
among the variables may be calculated by setting the left-hand side of 
equation ( 4 ) equal to the load on the water, letting £ = 0 , and proceeding 
as previously. It should also be noted that z = 0 for this case. Sub- 
stituting the above conditions leads to the following result: 
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z = C A,p 

b C v 2 G cp(A) sin t 


where 

C A ,p planing beam-loading coefficient, L^pgb 

L load on the water 

Cy speed coefficient, 



DISCUSSION 


The basic theory is believed to be adequate over the entire range 
of initial flight conditions, provided the proper two-dimensional water- 
mass variation is used. Since the water-mass variation of the present 
application is for deeply immersed chines, the application is limited to 
heavily loaded models having substantial chine immersion. Impacts with 
nonimmersed or moderately immersed chines have been analyzed in refer- 
ences 1 to 

In application of the present theory, it can be seen from equa- 
tions (21) and (23) that the value of G must be obtained in order 

to evaluate the coefficients from experimental data. This value depends 
upon the depth of immersion at the point of interest and the dead-rise 
angle. In the present application it is approximated by a constant. 

This approximation seems logical since the aspect-ratio correction for 
deep immersions asymptotically approaches 1. The value used to reduce 
the data of the present investigation was empirically determined for best 
fit; for the flat-bottom model having a range of C A from 36.5 to 133-7 

the value of G was found to be approximately 1. The data of 

references 4 and 5 obtained from tests at the Langley impact basin were 
used for comparison with the theory. These data covered a range of beam- 
loading coefficients C A from 36.5 to 133-7 tor the model with a flat 

bottom. 

The agreement between the theory and experiment for the accelera- 
tion coefficient is shown in figures 3 to 5* Figures 3 and 4 show the 
variation of the acceleration coefficient with the approach parameter at 
maximum load and maximum penetration, respectively, and figure 5 shows 
the variation of the acceleration coefficient with time. These figures 
show fairly good agreement between the experimental results and the theory. 
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The penetration coefficient at maxiimun penetration and maximum load 
is shown in figures 6 and respectively In general , the penetration 
is overestimated by theory, especially at the higher values of k. 

Considerable data for the body with a dead-rise angle of 30° and 
a range from 208 to 503 have been obtained at the Langley impact 

basin (ref. 4-). The trim angles ranged from 6° to 30° and the flight- 
path angles from approximately 2° to 20°. These data are compared with 
the present theory by using an empirically determined dead-rise factor 
to correct for the dead-rise effect. Thau is, the value for G cp(Ai) 

was taken as 0.6l instead of the 1 used for the case of 0° dead-rise 
angle. The results are shown in figure 8 for the maximum acceleration 
coefficient, and the agreement suggests that the empirically determined 
dead-rise constant is sufficient for correcting the data at the dead- 
rise angle of 30° and high values of beam-loading coefficient. As a 
first approximation, a linear variation of G cp( A-^ with dead-rise 

angle could be assumed to determine the intermediate values at these 
high beam-loading coefficients. However, more data are necessary to 
establish a function that will hold for all dead-rise angles. 

To see how far the present application of the theory could be 
extended to lower beam-loading coefficients, the data of references 4 
and 5 were used. These data were for flat-bottom models with C^ equal 

to 4.4 and 18.8, flight-path angles from L.79 to 21.21, and trim angles 
from 3° to 45° . It was found that fair agreement could be obtained for 
the maximum acceleration coefficient by again merely using a different 
constant for the value of G <p(Aq) with sach of the lower beam loadings . 

Figures 9 and 10 show the agreement for tie beam-loading coefficients 
of 18.8 and 4.4 when G q^Ap) is taken as 2 and 4, respectively. It 

is difficult to visualize any increase in <p(Ap) over the value of 1 
in any casej therefore, it appears that as C^ decreases to low values, 

G pg dependent on C ^ as well as on p. Of course, as C^ approaches 1 

the theory as presently applied would be expected to show errors because 
of the assumption of a two-dimensional water mass. 


CONCLUDING REMARKS 


A theory has been derived for the leads and motions of a deeply 
immersed prismatic body throughout a hydrodynamic impa-ct. The time and 
motion coefficients are presented in a generalized form, similar to that 
previously employed for the nonimmersed-chine case, which involves the 
use of an approach parameter k that depends only on the initial trim 
and flight-path angles. The use of this parameter reduces the number of 
independent variables and thereby simplifies presentation of results. 
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The theory is substantiated, over a wide range of initial flight 
conditions for bodies having a dead-rise angle of 0° and beam-loading 
coefficients C A that range from 36.5 to 133 -7- An empirically deter- 
mined factor was used in comparing the theoretical and experimental 
maximum impact loads for the body with a 30° dead-rise angle and beam- 
loading coefficients ranging from 208 to 530. The comparison showed 
good agreement. Fairly good agreement was also obtained for the maximum- 
load data of bodies with 0° dead-rise angle and beam-loading coeffi- 
cients of 18.8 and 4.4 by application of empirically determined correc- 
tion factors. 


Langley Research Center, 

National Aeronautics and Space Administration, 
Langley Field, Va., November 4, 1958. 
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Figure 3 .- Theoretical and experimental variation of maximum accelera- 
tion coefficient with approach parametei . p = 0°; b = 8 to 12 inches 

C A = 36-5 to 133 -7- 




Figure k.- Theoretical and experimental variation of acceleration coef 
ficient at maximum penetration with approach parameter. (3 = 0°; 
b = 8 to 12 inches: Ca = 36 . 5 to 133 -7* 
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Figure 8.- Theoretical and experimental variation of maximum accelera- 
tion coefficient with approach parameter for a body having a dead- 
rise angle of 30°. b = 5 inches; = 208 to 530; r = 6° to 30°; 

W = 933 to 2,350 pounds; 7 = 2.22° to 20.91; V Q = 3^-7 to 87.8 fps 
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Figure 10.- Theoretical and experimental variation of maximum accel 
eration coefficient with approach parameter for = 4 . 4 -. 0 = 1 

b = 20 inches; W = l,26l pounds. 
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